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Procede de traitement de donnees sonores et dispositif 
d 1 acquisition sonore mettant en oeuvre ce procede 

La presente invention concerne le traitement de donnees 
5 sonores . 

Des techniques relatives a la propagation d'une onde 
sonore dans l'espace tridiinensionnel r impliquant notamment 
une simulation et/ou une restitution sonores specialisees, 

10 mettent en oeuvre des procedes de traitement du signal 
audio appliques a la simulation de phenomenes acoustiques 
et psycho-acoustiques . De tels procedes de traitement 
pr^voient un encodage spatial du champ acoustique, .sa 
transmission et sa reproduction spatialisee sur pxi 

15 ensemble de hauts-parleurs ou sur des 6couteurs 
casque stereophonique . 

Parmi les techniques de son spatialise, on distingue deux 
categories de traitements complementaires l'une de l 1 autre 
20 mais qui sont g§neralement mise en oeuvre, l'une et 
1 1 autre , au sein d 1 un meme systdme 

D'une part, une premiere categorie de traitements concerne 
les procedes de synthese d'effet de salle, ou plus 

25 generalement d'effets environnementaux . A partir d'une 
description d'une ou plusieurs sources sonores (signal 
emis, position, orientation, direct ivite, ou autre) et en 
se basant sur un modele d'effet de salle (impliquant une 
geometrie de salle, ou encore une perception acoustique 

30 souhaitee) , on calcule et I'on decrit un ensemble de 
phenomenes acoustiques elementaires (ondes directes, 
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reflechies ou dif fractees) , ou encore un phenomene 
acoustique macroscopique (champ reverbere et dif f us) , 
permettant de traduire l'effet spatial au niveau d'un 
auditeur situe a un point choisi de perception auditive, 
dans l f espace tridimensionnel . On calcule alors un 
ensemble de signaux associes typiquement aux reflexions 
(sources " secondaires" , actives par re-£mission d'une onde 
principale recue, ayant un attribut de position spatiale) 
et/ou associes a une reverberation tardive (signaux 
decorreles pour un champ dif f us) . 

D 1 autre part, Une seconde categorie de procedes concerne 
le rendu positionnel ou directionnel de sources sonores. 
Ces precedes sont appliques a des signaux determines par 
un precede de la premiere categorie decrite ci-avant 
(impliquant des sources primaires et secondaires) en 
fonction de la description spatiale (position de la 
source) qui leur est associ<§e. En particulier, de tels 
procedes selon cette seconde categorie permettent 
d'obtenir des signaux k diffuser sur des hauts-parleurs ou 
ecouteurs, pour f inalement donner a un auditeur 
l 1 impression auditive de sources sonores placees a des 
positions respectives predetermines , autour de 
1« auditeur. Les proced6s selon cette seconde categorie 
sont qualifies de "createurs d r images sonores 
tridimensionnelles" , du fait de la repartition dans 
l ! espace tridimensionnel du ressenti de la position des 
sources par un auditeur. Des procedes selon la seconde 
categorie comportent generalement une premiere etape 
d'encodage spatial des evenements acoustiques elementaires 
qui produit une representation du champ sonore dans 
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l'espace tridimensional. Dans une seconde £tape, cette 
representation est transmise ou stockee pour un usage 
differe. Dans une troisieme etape, de decodage, les 
signaux decodes sont delivres sur des hauts-parleurs ou 
5 des ecouteurs d'un dispositif de restitution. 

La presente invention s'inscrit plutot dans la seconde 
categrorie precitee. Elle concerne en particulier 
l ! encodage spatial de sources sonores et une specification 

10 de la representation sonore tridimensionnelle de ces 
sources. Elle s 1 applique aussi bien a un encodage de 
sources sonores "virtuelles" (applications oil des sources 
sonores. sont simulees telles que des jeux, une conference 
spatialisee, ou autres) , qu'un encodage "acousfcigue" d'un ^ 

15 champ sonore naturel, lors d'une prise de son par un ou 

plusieurs reseaux tridimensionnels de microphones, i 

t 

Parmi les techniques envisageables de spatialisation du f 
son, 1 1 approche "ambisonique " est preferee . L 1 encodage 

20 ambisonique, qui sera decrit en detail plus loin, consiste 
a repr^senter des signaux relatifs a une ou plusieurs 
ondes sonores dans une base d 1 harmoniques spheriques (en 
coordonnees spheriques impliquant notamment un angle 
d 1 elevation et un angle azimutal, caracterisant une 

25 direction du ou des sons) . Les composantes representant 
ces signaux et exprimees dans cette base d' harmoniques 
spheriques sont aussi fonction, pour les ondes emises en 
champ proche, d'une distance entre la source sonore 
emettant ce champ et un point correspondant a l'origine de 

30 la base des harmonique spheriques. Plus particulierement , 
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cette dependance de la distance s'exprime en fonction de 
la frequence sonore, comme on le verra plus loin. 

Ce t te approche amb i sonique o f f re un grand nombr e de 
fonctionnalites possibles, notamment en terrae de 
simulation de sources virtuelles, et, de ntaniere generale, 
presente les avantages suivants : 

- elle traduit, de fagon rationnelle, la realite des 
phenomenes acoustiques et apporte un rendu auditif spatial 
realiste, convaincant et immersif ; 

- la representation des phenomenes acoustiques est 
scalable : elle off re une resolution spatiale qui peut 
etre adaptee & differentes situations. En effet, cette 
representation peut etre transmise et exploitee en 
fonction de contraintes de debit lors de la transmission 
des signaux encodes et/ou de limitations du dispositif de 
restitution ; 

- la representation ambisonique est flexible et il est 
possible simuler une rotation du champ sonore, ou encore, 
a la restitution, d 1 adapter le decodage des signaux 
ambisoniques & tout dispositif de restitution, de 
geometries diverses . 

Dans l 1 approche ambisonique connue, l'encodage des sources 
virtuelles est essentiellement directionnel . Les fonctions 
d'encodage reviennent a calculer des gains qui dependent 
de 1' incidence de l'onde sonore exprimee par les fonctions 
harmoniques sph^riques qui dependent de I 1 angle 
d' elevation et de 1' angle azimutal en coordonnees 
spheriques. En particulier, au decodage, on suppose que 
les hauts-parleurs, a la restitution, sont lointains . II 
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en resulte une distorsion (ou une incurvation) de la forme 
des fronts d'onde reconstruits. En effet, comme indique 
ci-avant, les composantes du signal sonore dans la base 
des harmoniques spheriques, pour un champ proche, 
5 dependent en fait aussi de la distance de la source et de 
la frequence sonore. Plus precisement, ces composantes 
peuvent s 1 exprimer mathematiquement sous la forme d l un 

pplynome dont la variable est inversement proportionnelle 

a la distance precitee et & la frequence sonore. Ainsi, 

10 les composantes ambisoniques , au sens de leur expression 
theorique, sont divergentes dans les basses frequences et, 
en particulier, tendent vers l'infini quand la frequence., 
sonore decroit vers zero, lorsqu'elles representent un son 
en champ proche emis par une source situee & une distance • Jl 

15 finie, Ce phenomene mathematique est connu, dans le L 
domaine . de la representation ambisonique, deja pour ,| 
l'ordre 1, par le terme de "bass boost", notamment .par : 
- M.A.GERZON, "General Metatheory of Auditory . 

Localisation", preprint 3306 of the 92 nd AES Convention, 

20 1992, page 52. 

Ce ph§nomene devient particulidrement critique pour des 
ordres d 1 harmoniques spheriques elevees impliquant des 
polynomes de puissance elevee . 

25 On connait par : 

SONTACCHI et HOLDRICH, "Further Inves tigations on 3D Sound 
Fields using- Distance Coding" (Proceedings of the COST G-6 
Conference on Digital Audio Effects (DAFX-01) , Limerick, 
Irlande, 6-8 Decembre 2001) , 



une technique pour prendre en compte une incurvation des 
fronts d f ondes au sein d'une representation proche d ! une 
representation ambisonique, dont le principe consiste a : 

- appliguer un encodage ambisonique (d'ordre eleve) aux 
signaux issus d'une prise de son virtuelle (simulee) , de 
type WFS (pour "Wave Field Synthesis") ; 

- et reconstruire le champ acoustique sur une zone d'aprds 
ses valeurs sur une frontiere de zone, se fondant ainsi 
sur le principe de HUYGENS - FRESNEL . 

Cependant, la technique presentee dans ce document, bien 
que prometteuse du fait qu'elle utilise une representation 
ambisonique a un ordre eleve, pose un certain nombre de 
probl^mes : 

- les ressources inf ormatiques necessaires pour le calcul 
de toutes les surfaces permettant d'appliquer le principe 
de HUYGENS - FRESNEL , ainsi que les temps de calcul 
necessaires, sont excessifs ; 

- des artefacts de traitement dits d 1 "aliasing spatial 11 
apparaissent a cause de la distance entre les microphones, 
a moins de choisir un maillage de microphone virtuels 
serre dans l f espace, ce qui alourdit les traitements ; 

- cette technique est dif f icilement transposable a un cas 
reel de capteurs i. disposer en reseau, en presence d'une 
source reelle, a 1 ' acquisition ; 

- a la restitution, la representation sonore 
tridimensionnelle est implicitement assujettie a un rayon 
fige du dispositif de restitution car le decodage 
ambisonique doit se faire, ici, sur un reseau de hauts- 
parleurs de memes dimensions que le reseau de microphones 
initial, ce document ne proposant aucun moyen d' adapter 
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l'encodage ou le decodage a d 1 autres tailles de 
dispositifs de restitution. 

Surtout, ce document presente un reseau horizontal de 
5 capteurs, ce qui suppose que les phenomenes acoustiques 
dont on tient compte, ici, ne se propagent que dans des 
directions, horizontales , ce qui exclut toute autre 
direction de propagation et qui, done, ne reprgsente pas 
la realite physique d f un champ acoustique ordinaire. 

10 

De . fagon .plus ggnerale, les techniques actuelles ne 
. permettent pas de traiter de fagon satisf aisante tout type 
de sources sonores , notamment en champ proche, mais plutot 
des sources sonores lointaines (ondes planes) , ce. qui 
15- correspond a une situation restrictive et artifioielle 
dans nombreuses applications. 

Un objet de la presente invention est de fournir un 
procede pour traiter, par encodage, transmission et 
20 restitution, tin type quelconque de champ sonore, en 
particulier 1 ' ef f et d'une source sonore en champ proche. 

Un autre objet de la presente invention est de fournir un 
procede permettant 1' encodage de sources virtuelles, non 
25 seulement en direction, mais aussi en distance, et de 
definir un decodage adaptable a un dispositif de 
restitution quelconque. 



30 



Un autre objet de la presente invention est de fournir un 
procede .de traitement robuste pour des sons de toutes 
frequences sonores (y compris les basses frequences) , 
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notamment pour la prise de son de champs acoustiques 
naturels a l'aide de reseaux tridimensionals de 
microphones . 

5 A cet effet, la presente invention propose un procede de 
traitement de donnees sonores, dans lequel : 

a) on code des signaux representatif s d'au moins un son se 
propageant dans 1 ' espace tridimensional et issu d'une 
source situee a une premiere distance d'un point de 

10 reference, pour obtenir une representation du son par des 
composantes exprimees dans une base d 1 harmoniques 
spheriques, d'origine correspondant audit point de 
reference, et 

b) on applique auxdites composantes une compensation d'un 
15 effet de champ proche par un filtrage qui est fonction 

d'une seconde distance definissant sensiblement , pour une 
restitution du son par un dispositif de restitution, une 
distance entre un point de restitution et un point de 
perception auditive. 

20 

Dans un premier mode de realisation, ladite source etant 
lointaine du point de reference, 

- on obtient des composantes d'ordres successifs m pour la 
representation du son dans ladite base d' harmoniques 

25 spheriques , et 

- on applique un filtre dont les coefficients, appliques 
chacun a une composante d ! ordre m, s ' expriment 
analytiquement sous la forme de 1 1 inverse d«un polynome de 
puissance m, dont la variable est inversement 

30 proportionnelle a la frequence sonore et a ladite seconde 



distance, pour compenser un effet de champ proche au 
niveau du dispositif de restitution. 

Dans un second mode de realisation, ladite source etant 
5 une source virtuelle prevue a ladite premiere distance, 

- on obtient des composantes d ! ordres successifs m pour la 
representation du son dans ladite base d ' harmoniques 
spheriques, et 

- on applique un filtre global dont les coefficients, 
10 appliques chacun a une composante d'ordre m, s ' expriment 

analytiquement sous la forme d'une fraction, dont : 

- le numerateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la ;. t 

frequence sonore et a ladite premiere distance, pour % 

;V 

15 simuler un effet de champ proche de la source a. 

virtuelle, et £ 

- le denominateur est un polynome de puissance m, £ 
dont la variable est inversement proportionnelle a la > 
frequence sonore et a ladite seconde distance, pour 

20 compenser l 1 effet du champ proche de la source 

virtuelle dans les basses frequences sonores . 

Pref erentiellement , on . transmet au, dispositif de 
restitution les dormees codees et filtrees aux etapes a) 
25 et b) avec un parametre representatif de ladite seconde 
distance - 

En complement ou en variante, le dispositif de restitution 
comportant des moyens de lecture d'un support memoire, on 
30 memorise sur un support memoire destine a etre lu par le 
dispositif de restitution les donnees codees et filtrees 
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aux etapes a) et b) avec un parametre representatif de 
ladite seconde distance. 

Avantageusement; prealablement a une restitution sonore 
par un dispositif de restitution comportant une pluralite 
de hauts-parleurs disposes a une troisieme distance dudit 
point de perception auditive, on applique aux donnees 
codees et filtrees un filtre d r adaptation dont les 
coefficients sont fonction desdites seconde et troisieme 
distances. 

Dans une realisation particuliere, les coefficients de ce 
filtre d' adaptation, appliques chacun a une composante 
d'ordre m, s 1 expriment analytiquement sous la forme d'une 
fraction, dont : 

- le numerateur est un polynome de puissance m, dont la 
variable est inversement proportionnelle a la frequence 
sonore et a ladite seconde distance, 

- et le denominateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la frequence 
sonore et a ladite troisieme distance. 

Avantageusement, pour la mise en oeuvre de l'etape b) , on 
prevoit : 

- pour des composantes d'ordre m pair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 
d'ordre deux ; et 

- pour des composantes d'ordre m impair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 
d'ordre deux et une cellule supplementaire d'ordre un. 
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Dans cette realisation, les coefficients d'un filtre 
audionumerique, pour une composante d'ordre m, sont 
definis a partir des valeurs numeriques des racines 
5 desdits polynomes de puissance m. 

Dans une realisation particuliere, les polynomes precites 
sont des polynomes de Bessel. 

10 A 1 ' acquisition des signaux sonores, on prevoit 
avantageusement un microphone comportant un r£seau de 
transducteurs acoustiques agences sensiblement sur la 
surface d'une sphere dont le centre correspond 

sensiblement audit point de reference, pour obtenir 

* I' 

15 lesdits signaux representatif s d'au moins un son :se 
propageant dans 1 1 espace tridimensional. ^ 

Dans cette realisation, on applique a l'etape b) un filtre 
global pour, d'une part, compenser un effet de champ 
20 proche en fonction de ladite seconde distance et, d' autre 
part, egaliser les signaux issus des transducteurs pour 
compenser une ponderation de directivite desdits 
transducteurs . 

25 Pref §rentiellement , on prevoit un nombre de transducteurs 
fonction d'un nombre total choisi de composantes pour 
representer le son dans ladite base d' harmoniques 
spheriques . 

30 Selon une caractiiristique avantageuse, on choisit 3. 
l'etape a) un nombre total de composantes dans la base des 
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harmoniques spheriques pour obtenir, £L la restitution, une 
region de l'espace autour du point de perception dans 
laquelle la restitution du son est fidele et dont les 
dimensions sont croissantes avec le nonabre total de 
■ 5 composantes . 

Pref erentiellement, on prevoit en outre un dispositif de 
restitution comportant un nombre de haut-parleurs au mo ins 
egal audit nombre total de composantes . 

10 

En variante, dans le cadre d'une restitution avec synthese 
binaurale ou transaurale : 

- on prevoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une 

15 distance choisie d'un auditeur, 

on obtient, pour cet auditeur, une information de 
ressenti attendu de la position dans l'espace de sources 
sonores situees a une distance de reference predeterminee 
de 1" auditeur pour 1 1 application d'une technique dite de 

20 "synthese binaurale" ou " transaurale " , et 

- on applique la compensation de l'etape b) avec ladite 
distance de reference sensiblement en tant que seconde 
distance . 

25 Dans une. variante ou l'on introduit une adaptation au 
dispositif de restitution a deux ecouteurs : 

- on prevoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une 
distance choisie d'un auditeur, 

30 - on obtient, pour cet auditeur, une information de 
ressenti de la position dans l'espace de sources sonores 
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situees a une distance de reference predeterminee de 
l'auditeur, et 

- prealablement a une restitution sonore par le dispositif 
de restitution, on applique aux donnees codees et filtrees 

5 aux etapes a) et b) un filtre d' adaptation dont les 
coefficients sont fonction de la seconde distance et 
sensiblement de la distance de reference. 

En particulier, dans le cadre d'une restitution avec 
10 synthese binaurale : 

- le dispositif de restitution comporte un casque a deux 
ecouteurs pour les oreilles respectives de l'auditeur, 

- et pref erentiellement , separement pour chaque ecouteur^ 
on applique le codage et le f iltrage des Stapes a) et bh 

15 pour des signaux respectifs destines a alimenter chaque- 
ecouteur, avec, en tant que premiere distance.** 
respectivement une distance separant chaque oreille d'une* 
position d ? une source §l restituer dans l'espace de 
restitution. ^ 

20 

Pref erentiellement , on met en forme, aux etapes a) et b) , 
un systeme matriciel comportant au moins : 

- une matrice comportant lesdites composantes dans la base 
des harmoniques spheriques, et 

25 - une matrice diagonale dont les coefficients 
correspondent a des coefficients de f iltrage de l r £tape 
b) , 

et on multiplie lesdites matrices pour obtenir une matrice 
resultat de composantes compensees . 

30 

De preference, a la restitution : 
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- le ciispositif de restitution comporte une pluralite de 
haut-parleurs disposes sensiblement a une meme distance du 
point de perception auditive, et 

- pour decoder lesdites donnees codees et filtrees aux 
5 etapes a) et b) et former des signaux adaptes pour 

alimenter lesdits haut-parleurs : 

* on forme un systeme matriciel comportant ladite 
_ matrice resultat de composantes compens§es et une 

matrice de d£codage pr^determinee, propre au 
10 dispositif de restitution, et 

* on obtient une matrice comportant des coefficients 
representatif s des signaux d 1 alimentation des hauts- 
parleurs par multiplication de la matrice resultat 
par ladite matrice de decodage* 

15 

La presente invention vise aussi un dispositif 
d' acquisition sonore, comportant un microphone muni d'un 
reseau de transducteurs acoustiques disposes sensiblement 
sur la surface d'une sphere. Selon 1' invention le 
20 dispositif comporte en outre une unite de traitement 
agencee pour : 

- recevoir des signaux 6manant chacun d'un transducteur , 

- appliquer auxdits signaux un codage pour obtenir une 
representation du son par des composantes exprimees dans 

25 une base d 1 harmonigues spheriques, d'origine correspondant 
au centre de ladite sphere, 

- et appliquer auxdites composantes un filtrage qui est 
fonction, d'une part, d'une distance correspondant au 
rayon de la sphere et, d 1 autre part, d'une distance de 

30 reference . 
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Pref erentiellement , le filtrage effectue par I 1 unite de 
traitement consiste, d'une part, a §galiser, en fonction 
du rayon de la sphere, les signaux issus des transducteurs 
pour compenser une ponderation de directivite desdits 
transducteurs et, d* autre part, a compenser un effet de 
champ proche en fonction de ladite distance de reference. 

D^autres avantages et caracteristiques de 1» invention 
apparaitront a la lecture de la description detaillee ci~ 
apres et a l'examen des figures qui 1 1 accompagnent , sur 
lesquelles : 

-* la figure 1 illustre schematiquement un syst^me 
d 1 acquisition et creation, par simulation de sources, 
virtuelles, de signaux sonores, avec encodage, 
transmission, decodage et restitution par un dispositif de 
restitution spatialise, 

- la figure 2 represente plus precisement un encodage de 
signaux definis k la fois en intensite et par rapport a la 
position d'une source dont ils sont issus, 

- la figure 3 illustre les parametres en jeu dans la 
representation ambisonique, en coordonnees spheriques ; 

- la figure 4 illustre une representation par une metrique 
tridimensionnelle dans un rep^re de coordonnees 

spheriques, d 1 harmoniques spheriques Y° m de differents 
ordres ; 

- la figure 5 est un diagramme des variations du module de 
fonctions radiales j m (kr) , qui sont des fonctions de 
Bessel spheriques, pour des valeurs d'ordre m successives, 
ces fonctions radiales intervenant dans la representation 
ambisonique d»un champ de pression acoustique ; 
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- la figure 6 represents 1 1 amplification due a l'effet de 
champ proche pour differents ordres successifs m, en 
particulier dans les basses frequences ; 

- la figure 7 represente schematiquement un dispositif de 
5 restitution comportant une pluralite de hauts-parleurs 

HPi , avec le point (reference P) de perception auditive 
precite, la premiere distance pr^citee (referencee p) et 
la - seconde distance, pre.citee . (referencee R> 

- la figure 8 represente schematiquement les parametres 
10 mis en jeu dans l f encodage ambisonique, avec un encodage 

directionnel, ainsi qu'un encodage de distance selon 
1' invent ion ; 

- la figure 9 represente des spectres d'energie des 
filtres de compensation et de champ proche simules pour 

15 une premiere distance d'une source virtuelle p = 1 m et 
une pr<S-compensation de hauts-parleurs situes a une 
seconde distance R = 1,5 m ; 

- la figure 10 represente des spectres d 1 energie des 
filtres de compensation et de champ proche simules pour 

20 une premiere distance de la source virtuelle p = 3 m et 
une pre-compensation de hauts-parleurs situes a une 
distance R = 1, 5 m ; 

- la figure 11A represente une reconstruction du champ 
proche avec compensation, au sens de la presente 

25 invention, pour une onde spherique dans le plan 
horizontal ; 

- la figure 11B, a comparer avec la figure 11A, represente 
le front d'onde initial, issu d*une source S ; 

- la figure 12 represente schematiquement un module de 
30 filtrage pour adapter les composantes ambisoniques regues 

et pre-compensees a l 1 encodage pour une distance de 
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reference R en tant que seconde distance, a un dispositif 
de restitution comportant une pluralite de hauts-parleurs 
disposes a une troisieme distance R 2 d'un point de 
perception auditive ; 

- la figure 13A represente schematiquement la disposition 
d ! une source sonore M f a la restitution, pour un auditeur 
utilisant un dispositif de restitution appliquant Tine 
synthese binaurale, avec une source erne t tant en champ 
proche ; 

- la figure 13B represente schematiquement les etapes 
d'encodage et de decodage avec effet de champ proche dans 
le cadre de la synthese binaurale de la figure 13 A a 
laquelle est combine un encodage/decodage ambisonique ; 

- la figure 14 represente schematiquement le traitement 
des signaux issus d'un microphone comportant une plurality 
de capteurs de pression agenc£s sur une sphere, 3. titre 
illustratif, par encodage ambisonique, egalisation et 
compensation de champ proche au sens de l 1 invention. 

On se refere tout d f abord a la figure 1 qui represente a 
titre illustratif un systeme global de spatialisation 
sonore. Un module la de simulation d'une scdne virtuelle 
definit un objet sonore comme une source virtuelle d'un 
signal, par exemple mohophonique , de position choisie dans 
1 1 espace tridimensional et qui definit une direction du 
son. II peut etre prevu en outre des specifications de la 
geometrie d'une salle virtuelle, pour simuler une 
reverberation du son. Un module de traitement 11 applique 
une gestion d'une ou plusieurs de ces sources par rapport 
a un auditeur (definition d'une position virtuelle des 
sources par rapport a cet auditeur) . II met en ceuvre un 
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processeur d'effet de salle pour simuler des 
reverberations ou autres en appliguant des retards et/ou 
des filtrages usuels . Les signaux ainsi construits sont 
transmis a un module 2a d'encodage spatial des 
5 contributions elementaires des sources. 

Parallelement, une prise de son naturelle peut etre 
effectuee dans le cadre d'un enregistrement sonore par un 
ou plusieurs microphones disposes de fagon choisie par 

10 rapport aux sources reelles (module lb) . Les signaux 
captes par les microphones sont encodes par un module 2b. 
Les signaux acquis et encodes peuvent etre transformes 
selon un format de representation intermediaire (module 
3b), avant d'etre mixes par le module 3 aux signaux 

15 generes par le module la et encodes par le module 2a (issu 
des sources virtuelles) . Les signaux mixes sont ensuite 
transmis, ou encore memorises sur un support, en vue d f une 
restitution ulterieure (fleche TR) . lis sont ensuite 
appliques a un module de decodage 5, en vue de la 

20 " restitution sur un dispositif de restitution 6 comportant 
des hauts-parleurs. Le cas echeant, 1 1 etape de decodage 5 
peut §tre precedee d 1 une etape de manipulation du champ 
sonore, par exemple par rotation, grace i. un module de 
traitement 4 prevu en amont du module de decodage 5 . 

25 

Le dispositif de restitution peut se presenter sous la 
forme d'une multiplicity de hauts-parleurs, agenc^s par 
exemple a la surface d J une sphere dans une configuration 
tridimensionnelle (periphonique) pour assurer, a la 
30 restitution, notamment un ressenti d'une direction du son 
dans 1 1 espace tridimensional . A cet effet, un auditeur se 
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place generalement au centre de la sphere formee par le 
reseau de haut-parleurs, ce centre correspondant au point 
de perception auditive cite ci-avant. En variante, les 
hauts-parleurs du dispositif de restitution peuvent etre 
agences dans un plan (configuration panoramique 
bidimensionnelle) , les hauts-parleurs etant disposes en 
particulier sur un cercle et 1' auditeur se plagant 
habituellement au centre de ce cercle. Dans une autre 
variante, le dispositif de restitution peut se presenter 
sous la forme d'un dispositif de type "surround" (5.1). 
Enfin, dans une variante avantageuse, le dispositif de 
restitution peut se presenter sous la forme d'un casque a 
deux ecouteurs pour une synthese binaurale du ^ :son 
restitue, qui permet a 1> auditeur de ressentir ,,^ne 
direction* des sources dans l'espace tridimensional, .cope 
on le verra . plus loin de fa(?oii d6taillee. Un jtel 
dispositif de restitution a deux hauts-parleurs, pour;, un 
ressenti dans l'espace tridimensionnel, peut se presenter 
aussi sous la., forme d'un dispositif de restitution 
transaurale, a deux hauts-parleurs disposes a une distance 
choisie d'un auditeur. 

On se refere maintenant a la figure 2 pour decrire un 
encodage spatial et un decodage pour une restitution 
sonore tridimensionnelle, de sources sonores §lementaires . 
On transmet a un module d l encodage spatial 2 le signal 
issu d'une source 1 a N, ainsi que sa position (reelle ou 
virtuelle) . Sa position peut etre aussi bien definie en 
terme d' incidence (direction de la source vue de 
i' auditeur) qu'en terme de distance entre cette source et 
un auditeur. La plural ite des signaux ainsi encodes permet 
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d'obtenir une representation multi-canale d'un champ 
sonore global. Les signaux encodes sont transmis (fleche 
TR) a un dispositif de restitution sonore 6, pour une 
restitution sonore dans I'espace tridimensionnel, comme 
5 indique ci-avant en reference a la figure 1. 

On se refere maintenant a la figure 3 pour decrire ci- 

apres la representation ambisonique par des harmoniques 

spheriques dans 1 1 espace tridimensionnel, d'un champ 
10 acoustique. On considere une zone autour d'une origine O 
(sphere de rayon R) exempte de source acoustique. On 
adopte un systeme de coordonnees spheriques dans lequel 
chaque vecteur 7 des 1 1 origine O a un point de la sphere 
est deer it par un azimut 9 r , une elevation 5 r et un rayon r 
15 (correspondant a la distance a I'origine O) . 

Le champ de pression p{7) a 1 1 interieur de cette sphere 
(r < R ou R est le rayon de la sphere) peut s 1 ecrire dans 
le domaine frequentiel comme une serie dont les termes 

20 . sont les produits ponderes de fonctions angulaires y„ in {d,S) 
et de fonction radiale j m (kr) qui dependent ainsi d'un 
terme de propagation ou k=27tf/c / ou f est la frequence 
sonore et c est la vitesse du son dans le milieu de 
propagation . 

25 

Le champ de pression s 1 exprime alors par : 

p{r)=f j f'j m (kr) 2 ^ nn Y^ Nm {O r ,S r ) [Al] 
m=0 0<n<>m,cr=±} 
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L 1 ensemble des facteurs de ponderation B? nn , qui sont 
implicitement fonction de la frequence, decrivent ainsi le 
champ de pression dans la zone consideree. Pour cette 
raison, ces facteurs sont appeles "composemtes harmoniques 
spheriques" et representent une expression frequentielle 
du son (ou du champ de pression) dans la base des 

harmoniques spheriques Y% n - 

Les fonctions angulaires sont appelees "harmoniques 
sphSriques " et sont definies. par : 



ou 

P mn (s±nS) sont des fonctions de Legendre de degr6 m et 
d ■ ordre n ; 

5 P/ g est le symbole de Kronecker (egal 1 si p=q et 0, 
sinon) 

Les harmoniques spheriques forment une base orthonormee ou 
les produits scalaires entre composantes harmoniques et, 
de facjon generale entre deux fonctions F et G, sont 
respect ivement definies par : 



[A2J 



V \m + n 



cos n6 si a = +1 
sin nO si a = —l 




[A'2] 



(F\G) 4/r =^F(^)G(^),W5) 

Les harmoniques spheriques sont des fonctions reelles 
borneeS; comme represents sur la figure 4, en fonction de 
l 1 ordre m et des indices n et a. Les parties sombres et 
claires correspondent respectivement aux valeurs positives 
et negatives des fonctions harmoniques spheriques. Plus 
l 1 ordre m est eleve et plus la frequence angulaire (et 
done la discrimination entre fonctions) est elevee. Les 
fonctions radiales j m (kr) sont des fonctions de Bessel 
spheriques, dont le module est illustre pour quelques 
valeurs de l 1 ordre m dans la figure 5. 

On peut donner une interpretation de la representation 
ambisonique par une base d 1 harmoniques spheriques comme 
suit. Les composantes ambisoniques de metne ordre m 
expriment finalement des "derivees" ou des "moments" 
d 1 ordre m du champ de pression au voisinage de l'origine O 
(centre de la sphdre representee sur la figure 3) . 

En particulier, B$ = W decrit la grandeur scalaire de la 

pression, tandis que B\\ = X y fifj 1 = Y y B+o = Z sont lies aux 
gradients de pression (ou encore a la velocite 
particulaire) a l'origine O. Ces quatre premieres 
composantes W, X, Y et Z sont obtenues lors d'une prise de 
son naturelle a l'aide de microphones omnidirectif s (pour 
la composante W d 1 ordre 0) et bidirectifs (pour les trois 
autres composantes suivantes) . En utilisant un plus grand 
nombre de transducteurs acoustiques, un traitement 
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approprie, notamment par egalisation, permet d ! obtenir 
d ! avantage de composantes ambisoniques (ordres m plus 
eleves superieurs a 1) . 

5 En prenant en compte des composantes supplementaires 
d ! ordre plus eleve (superieur a 1), done en augment ant la 
resolution angulaire de la description ambisonique, on 
accede a „H ne approximation du champ de pression sur un 
voisinage plus large au regard de la longueur d'onde de 

10 l ! onde sonore, autour de 1 1 origine O. On comprendra ainsi 
qu'il existe une relation etroite entre la resolution 
angulaire (ordre des harmoniques spheriques) et la portee 
radiale (rayon r) qui peut itre representee. En bref 
lorsque l'on s 1 ecarte spatialement du point d 1 origine O de , 

15 la figure 3, plus le nombre de composantes ambisoniques, § 
est eleve (ordre M eleve) et meilleure .. est la4 
representation du son par 1 1 ensemble de ces composantes. « 
• ambisoniques. On comprendra aussi que la representation: 
ambisonique du son est toutefois moins satisf aisante au 

20 fur et a mesure que l'on s 1 eloigne de 1' origine O. Get 
effet devient critique en particulier pour des frequences 
sonores elevees (de longueur d f onde courte) . On a done 
interet a obtenir un nombre de composantes ambisoniques 
qui soit le plus grand possible, ce qui permet de creer 

25 une region de 1 1 espace autour du point de perception, dans 
laquelle la restitution du son est fidele et dont les 
dimensions sont croissantes avec le nombre total de 
composantes . 

30 On decrit ci-apres une application a systeme 
d 1 encodage/ transmission/restitution d'un son spatialise. 




24 



En pratique, un systeme ambisonique prend en compte un 
sous- ensemble de composantes harmoniques spheriques, comme 
decrit ci-avant. On parle d'un systeme d'ordre M lorsque 

5 celui-ci prend en compte des composantes ambisoniques 
d'indice m < M. Lorsqu f il s'agit d'une restitution par un 
dispositif de restitution a hauts-parleurs, on comprendra 
que si ces hauts-parleurs sont disposes dans un plan 
horizontal, seules les harmoniques d f indice m=n sont 

10 exploitees. En revanche, lorsque le dispositif de 
restitution comporte des hauts-parleurs disposes sur la 
surface d'une sphere ("periphonle" ) , on peut en principe 
exploiter autant d'harmoniques qu'il existe de haut- 
parleurs . 

15 

On designe par la reference S le signal de pression porte 
par une onde plane et capte au point O correspondant au 
centre de la sphere de la figure 3 (origine de la base en 
coordormees spheriques) . L 1 incidence de l f onde est decrite 
20 par l'azimut 9 et l 1 elevation 5. L 1 expression des 
composantes du champ associe a cette onde plane est donnee 
par la relation : 



B? nn =SY? m WJ) [A3] 

25 

Pour encoder (simuler) une source en champ proche a une 

distance p de 1» origine O, on applique un filtre F^ pfc) 
pour "incurver" la forme des fronts d f onde, en considerant 
qu'un champ proche emet, en premiere approximation, une 
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onde • spherique. Les composantes encodees du champ 
deviennent : 



et 1' expression du filtre precite Fjf lc ^ est donnee par la 
relation : 



10 ou a> = 27cf est la pulsation de l'onde, f etant la 
frequence du son. 

Ces deux dernieres relations [A4] et [A5] montrent & 
finalement que, aussi bien pour une source virtuelle % 

15 (simul<§e) que pour une source reelle en champ proche, les >- 
composantes du son dans la representation ambisonique ..- 
s'expriment mathematiquement (en particulief r 

analytiquement ) sous la forme d'un . polyndme, ici de 
Bessel; de puissance m et dont la variable (c/2jcop) est 

20 inversement proportionnelle a la frequence sonore. 

Ainsi, on comprendra que : 

- dans le cas d'une onde plane, 1 r encodage produit des 
signaux qui ne different du signal d'origine que d'un gain 

25 reel, fini, ce qui correspond a un encodage purement 
directionnel (relation [A3] ) ; 

- dans le cas d'une onde spherique (source en champ 
proche) , le filtre supplementaire F^ p/c \co) encode 
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1 1 information de distance en introduisant , dans 
1' expression des composantes ambisoniques, des rapports 
d 1 amplitudes complexes qui dependent de la frequence, 
comme exprime dans la relation [A5] , 

5 

II est a noter que ce filtre supplementaire est de type 
"integrrafceur avec un effet d ! amplification croissant et 
divergent (non-borne) au fur et a mesure que les 
frequences sonores decroissent vers zero. La figure 6 

10 montre, pour chaque ordre m, une augmentation du gain en 
basses frequences (ici la premiere distance p = lm) . II 
s'agit done de filtres instables et divergents lorsque 
l'on cherche a les appliquer a des signaux audio 
quelconques. Cette divergence est d'autant plus critique 

15 pour les ordres m de valeur elevee. 

On comprendra en particulier, a partir des relations [A3] , 
[A4] , et [A5] , que la modelisation d'une source virtuelle 
en champ proche presente des composantes ambisoniques 

20 divergentes en basses frequences, de fagon 

particulierement critique pour des ordres m eleves, comme 
represents sur la figure 6. Cette divergence, dans les 
basses frequences, correspond au phenomene de "bass boost" 
enonce ci-avant. II se manifeste aussi en acquisition 

25 sonore, pour des sources reelles . 

Pour cette raison notamment, 1 1 approche ambisonique, en 
particulier pour des ordres m eleves, n'a pas connu, dans 
l'etat de la technique, une application concrete (autre 
30 que theorique) dans le traitement du son. 
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On comprend en pax-ticulier qu'une compensation du champ 
proche est necessaire pour respecter, a la restitution, la 
forme des fronts d'ondes encodes dans la representation 
ambisonique. En se referant a la figure 7, un dispositif 
5 de restitution comporte une pluralite de hauts-parleurs 
HP ±/ disposes a une meme distance R, dans l'exemple 
decrit, d'un point de perception auditive P. Sur cette 
figure 7 : 

chaque point ou se situe un haut-parleur HP ± correspond 
10 a un point de restitution enonce ci-avant, 

- le point P est le point de perception auditive enonce 
ci-avant , 

- ces points sont separes de la seconde distance R enoncee , ? 
ci-avant; 

■ 

15 tandis que sur la figure 3 decrite ci-avant : 

- le point O correspond au point de reference, enonce ci- l v 
avant, qui forme l'origine de la base des harmoniques 
spheriques, 

- le point M correspond & la position d'une source (reelle 
20 ou virtuelle) situee a la premiere distance p, enoncee ci- 

avant, du point de reference O. 

Selon 1' invention, on introduit une pre -compensation du 
champ proche au stade meme de l ! encodage, cette 
25 compensation mettant en jeu des filtres de la forme 

analytique 1 et qui s'appliquent aux composantes 

F m \fi>) 

ambisoniques B° ln precitees . 



. — . 1 
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Selon l f un des avantages que procure 1 ' invention, 

1 'amplification Fl p/c \co) dont l f effet apparait sur la 
figure 6 est compensee par 1 1 attenuation du filtre 

applique des l f encodage — — ^ • En particulier, les 

5 " coefficients de ce filtre de compensation — rsrr^ sont 

croissants avec la frequence du son *'et,' en* 'particulier ," 
tendent . vers zero, pour les basses frequences. 
Avantageusement , cette pre -compensation, effectuee des 
l'encodage, assure que les donnees transmises ne sont pas 
10 divergentes pour les basses frequences. 

Pour indiquer la signification physique de la distance R 
qui intervient dans le filtre de compensation, on 
considere, a titre illustratif , une onde plane r^elle, 

15 initiale, k 1 1 acquisition des signaux sonores . Pour 
simuler un effet de champ proche de cette source 
lointaine, on applique le premier filtre de la relation 
[A5] , comme indique dans la relation [A4] . La distance p 
represente alors une distance entre une source virtuelle 

20 proche M et le point O representant 1 1 origine de la base 
spherique de la figure 3. On applique ainsi un premier 
filtre de simulation de champ proche pour simuler la 
presence d'une source virtuelle a la distance p decrite 
ci-avant. Neanmoins, d'une part, comme indique ci-avant, 

25 les termes du coefficient de ce filtre divergent dans les 
basses frequences (figure 6) et, d 1 autre part, la distance 
p precitee ne representera pas forcement la distance entre 
les hauts-parleurs d*un dispositif de restitution et un 
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point P de perception (figure 7), Selon 1" invention, on 
applique une pre-compensation, a l'encodage, mettant en 



qui permet, d'une part, de transmettre des signaux bornes, 
et, d r autre part, de choisir la distance R, des 
l'encodage, pour la restitution du son a partir des hauts- 
parleurs HPi, tel que, represents sur la figure 7. En 
particulier, on comprendra que si l'on a simule, 3. 
1 'acquisition, une source virtuelle plac<§e a la distance p 
de l'origine 0, §. la restitution (figure 7), un auditeur 
place au point P de perception auditive (& une distance R 
des hauts -parleurs HP ± ) ressentira, a l'audition, la 
presence d'une source sonore S, placSe a la distance p du 
point de perception P et qui correspond a la source 
virtuelle simulee lors de 1 ' acquisition. 

Ainsi, la pre -compensation du champ proche des hauts- 
parleurs (places a la distance R) , au stade de l'encodage, 
peut etre combinee a un effet de champ proche simule d'une 
source virtuelle placee a une distance p. A l'encodage, on 
met finalement en jeu un filtre total resultant, d T une 
part, de la simulation . du champ proche, et, d' autre part, 
de la compensation du champ proche, les coefficients de ce 
filtre pouvant s'exprimer analytiquement par la relation : 



jeu un filtre de type 



comme indique ci-avant, ce 



H NFC(p/c.R,c) (co) = 



F { m RIC \o» 



[All] 
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Le filtre total donne par la relation [All] est stable et 
constitue la partie "encodage de distance" dans 1 1 encodage 
ambisonique spatial selon 1' invention, tel que represents 
sur la figure 8. Les coefficients de ces filtres 
5 correspondent a des f one t ions de trans fert monotones de la 
frequence, qui tendent vers la valeur 1 en hautes 
frequences et vers la valeur (R/p) m en basses frequences. 
. ..En .se .referant a „ la., figure .9, les spectres^ d'£nergie des 

filtres H% FC{plc > Rtc) ((D) traduisent 1 1 amplification des 
10 composantes encodees, dues a l'effet de champ de la source 
virtuelle (placee ici a une distance p = 1 m) , avec une 
pr§- compensation du champ des hauts-parleurs (places a une 
distance R = 1,5 m). L 1 amplification en decibels est done 
positive lorsque p < R (cas de la figure 9) et negative 
15 quand p > R (cas de la figure 10 ou p = 3 m et R = 1,5m). 
Dans un dispositif de restitution spatialisee, la distance 
R entre un point de perception auditive et les haut- 
parleurs HPi est ef f ectivement de I'ordre de un ou 
quel que s metres, 

20 

En se refer ant a nouveau a la figure 8, on comprendra que, 
outre les parametres de direction 6 et 8 habituels, on 
transmettra une information sur les distances qui 
interviennent a 1' encodage. Ainsi, les fonctions 

25 angulaires correspondant aux harmoniques spheriques Y^ 0%S ^ 
sont conservees pour 1 1 encodage directionnel . 

Toutefois, au sens de la presente invention, on prevoit en 
outre des filtres totaux (compensation de champ proche et, 
30 le cas echeant, simulation d'un champ proche) 
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Hm FCi<p/CiR/c) (co) qui sont appliques aux composantes 
ambisoniques, en fonction de leur orcire m, pour realiser 
1' encodage de la distance, comme represents sur la figure 
8 . Un mode de realisation de ces f iltres dans le domaine 
audionumerique sera decrit en detail plus loin. 



On remarquera en particulier que ces filtres peuvent etre 
appliques des meme 1 * encodage de distance (r) et avant 
■meme I 1 encodage de direction (9,8). On comprendra ainsi 
que les etapes a) et b) ci- avant peuvent etre rassemblees 
en une meme etape globale, ou meme etre interverties (avec 
un encodage de distance et filtrage de compensation,., 
suivis d'un encodage de direction) . Le procede selon- : 
I 1 invention ne se limite done pas k une mise en oeuvre :tj 
successive dans le temps des etapes a) et b) . 

La figure 11A represente une visualisation (vue de dessus)* 
d'une reconstruction d'un champ proche avec' compensation, 
d'une onde spherique, dans le plan horizontal (avec les 
memes parametres de distance que ceux de la figure 9) , 
pour un systeme d'ordre total M = 15 et une restitution 
sur 32 hauts-parleurs . Sur la figure 11B, on a represente 
la propagation de 1 1 on,de sonore initiale a partir d'une 
source en champ proche situee & une distance p d'un point 
de l'espace d 1 acquisition qui correspond, dans 1 ' espace de 
restitution, au point P de la figure 7 de perception 
auditive. On remarque sur la figure 11A que les auditeurs 
(symbolises par des tetes schernatisees) peuvent localiser 
la source virtuelle en un meme lieu geographique situe a 
la distance p du point de perception P sur- la figure 11B. 
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On verifie bien ainsi que la forme du front d'onde encode 
est respectee apres decodage et restitution. Toutefois, on 
constate sensiblement des interferences a droite du point 
P tel que represents sur la figure 11A qui sont dues au 
fait que le nombre de hauts-parleurs (done de composantes 
ambisoniqTj.es prises en compte) n'est pas suffisant pour 
restituer parfaitement le front d'ondes en jeu sur toute 
la surface delimitee par les haut-parleurs. 

Dans ce qui suit, on decrit, a titre d'exemple, 
l'obtention d ! un filtre audionumerique pour la mise en 
ceuvre du procede au sens de l f invention. 

Comme indique ci-avaxit, si l'on cherche a simuler un effet 
de champ proche, compense des l'encodage, on applique aux 
composantes ambisoniques du son un filtre de la forme : 



H NFC {P lc,R, c) } . F >n P,C) ^ [All] 

De 1' expression de la simulation d'un champ proche donne 
par la relation [A5] , il apparait que pour des sources 
lointaines (p = co) , la relation [All] devient simplement : 

^ H NFCM/c) (0)) [A12] 



F<« /c hco) m 



II apparait done de cette derniere relation [A12] que le 
cas ou la source a simuler emet en champ lointain (source 
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lointaine) n'est qu ! un cas particulier de 1 1 expression 
generale du filtre formulee dans la relation [All] . 

Dans le domaine des traitements audi onumeri que s , un 
procedeS avantageux pour definir un filtre numerique a 
partir de I 1 expression analytique de ce filtre dans le 
domaine analogique a temps continu consiste en une 
" trans formee bilineaire " . 

On exprime d'abord la relation [A5] sous la forme d'une 
transformee de Laplace, ce qui correspond a : 



ou x = p/c (c etant la vitesse acoustique dans le milieu, 
typiquement 340 m/s dans l'air) . 

La transformee bilineaire consiste a presenter, pour une 
frequence d 1 echantillonnage f s , la relation [All] sous la 
forme : 

si m est impair et 




[A13] 




34 



si m est pair, 



1- 



-1 



ou z est defini par p = 2f s par rapport a la relation 

[A13] precedente, 



et avec 



a 



~ _2 



1 - 



X 



a 



2 " 



et 



X~2 = A + ~ 1 ~ 



10 



*0 



a 



et *> +1 > /2 =- 



1 + 



X 



m,q 



a 



15 



ou a = 4 f s R/ c pour x=a 
et a = 4f s p/c pour x=b 

X ra , q sont les g racines successives du polyn6me de Bessel 



p - V ("* + ») ? Y m-n 



-n(' v --v») 
</=i 
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et sont exprimees dans le tableau 1 ci-apres, pour 
differents ordres m, sous les formes respectives de leur 
partie reelle, leur module (separes par une virgule) et 
leur valeur (reelle) lorsque m est impair. 
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Tableau 1 : valeurs R t \X m ^ q '\, X mq \ (et R t .[X mm ) lorsque m 
est impair) d'un polynSme de Bessel calculees a l'aide du 
logiciel de calcul MATLAB© . 



m=l 


-2 . 0000000000 


m=2 


-3.0000000000, 3.4641016151 


m=3 


-3.6778146454, 5.0830828022 ; -4.6443707093 


A 

m=4 


a Ofi7^7ft*7Qdd £ 77fl7^1 Rft^4 • ^ 7Q7471 90^fi Q46R29877S 

1 , f J / O / 7"^ / D. r 'O / J X J O J 1 , Z> . / j^liii^v JO | O . V *± O .<£ Zf O / / 


m=5 


-4.6493486064, 8.5220456027 ; -6.7039127983, 7.5557873219 ; 
-7.2934771907 


m=6 


-5.0318644956, 10.2983543043 ; -7.4714167127, 9.1329783045 ; 
-8 .4967187917, 8 .6720541026 


m=7 


-5.3713537579, 12.0990553610 ; -8.1402783273, 10.7585400670 ; 
-9.5165810563, 10.1324122997 ; -9.9435737171 


m=8 


-5.6779678978, 13.9186233016 ; -8.7365784344, 12.4208298072 ; 
-10.4096815813, 11.6507064310 ; -11.1757720865, 11.3096817388 


m=9 


-5.9585215964, 15.7532774523 ; -9.2768797744, 14.1121936859 ; 
-11.2088436390, 13.2131216226 ; -12.2587358086, 12.74194143 92 ; 
-12 .5940383634 


m=10 


-6.2178324673, 17.6003068759 ; -9.77243 91337, 15.8272658299 ; 
-11.9350566572, 14.8106929213 ; -13.2305819310, 14.2242555605 ; 
-13.8440898109, 13 . 9524261065 


m=ll 


-6.4594441798, 19.4576958063 ; -10.2312965678, 17.5621095176 ; 
-12.6026749098, 16.4371594 915 ; -14.1157847751, 15 ..7463731900 ; 
-14.9684597220, 15.3663558234 ; -15.2446796908 


m=12 


-6.6860466156, 21.323 9012076 ; -10.6594171817, 19.3137363168 ; 
-13.2220085001, 18.0879209819 ; -14.9311424804, 17.3 012295772 ; 
-15.9945411996, 16.824216 5032 ; -16.5068440226, 16.5978151615 


m=13 


-6.8997344413, 23.1977134 580 ; -11.0613619668, 21.0798161546 ; 
-13.8007456514, 19.7594692366 ; -15.6887605582, 18.88367673 5 9 ; 
-16.9411835315, 18.3181073 534 ; -17.6605041890, 17.9988179873 ; 
-17 .8554193236 


m=l4 


-7.1021737668, 25.0781652657 ; -11.4407047669, 22.8584 924996 ; 
-14.3447919297, 21.4490520815 ; -16.3976939224, 20.4898067617 ; 
-17.8220011429, 19.8423306934 ; -18.7262916698, 19.4 38 9130000 ; 
-19.1663428016, 19.244 7495545 
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m-15 


-7.2947137247, 26.9644699653 ; 
-14.8587939669, 23.1544615283 
-18.6471986915, 21.3925954403 
-20.3418287818, 20.6361378957 


-11.8003034312, 24.6482592959 ; 
; -17.0649181370, 22.1165594535 ; 
• -19.7191341042, 20.9118275261 ; 
? -20.5462183256 


m-16 


-7.4784635949, 28.8 5597844 87 
-15.3464816324, 24 . 873 8935490 
-19.4246523327, 22 . 9655586516 
-21.4379698156, 22.0627133 056 


-12.1424827551, 26.447876 0 957 ; 

• -17.6959363478, 23.76147 99683 ; 

• -20.6502404436, 22.4128776078 ; 

• -21 .8237730778, 21 . 89266 62470* 


m=17 


-7.6543475694, 30.7521483222 ; 
-15. 81089906 91, 26 . 605 8519104 
-20Yl605894^729, '2~4 . 5535534450 
-22.4668764601, 23 . 5193877036 
-23.1970582109 


-12.4691619784, 28.2563077987 ; 

• -18.2951775164, 25.4225585034 
? —21.5282660840,'23.9384287933 

• -23 . 0161527444, 23 . 2766166711 




m=18 


-7.8231445835, 32.6525213363 ; 
-16. 2545681590, 28.34 90792784 
-20. 8600257104 , 26 . 1693913642 
-23 .4378933084, 25.0022244227 
-24 .4798038436, 24 . 5412441597 


-12,7819455282, 30.0726807554 ; 

• -18.8662638563, 27.0981271991 

• -22.3600808236, 25.4856138632 

• -24.1362741870, 24.6925542646 




m=19 


-7.9855178345, .34.5567065132 ; 
-16.6796008200, 30.1025072510 , 
-21.5270719955, 27.7962699865 , 
-24.3584393996, 26,5081174988 , 
-25.6855663388, 25.9191817486 ( 


-13.0821901901, 31.8962504142 ; 

-19.4122071436, 28.786 7778706 , 
• -23.1512112785, 27.0520753105 

-25.1941793616, 26.13 63057951 

-25.8480312755 





On implemente ainsi les f litres nuriieriques , a partir des 
valeurs du tableau 1, en prevoyant des cascades de 
cellules d'ordre 2 (pour m pair), et une cellule 
supplemental re (pour m impair) , a partir des relations 
[A14] donnees ci-avant . 

On realise ainsi des filtres numeriques sous une forme de 
reponse impulsionnelle infinie, aisement parametrable 
comme montre ci-avant. II est a noter qu'une 
implementation sous une forme de reponse impulsionnelle 
f inie peut §tre envisagee et . consiste a calculer le 
spectre complexe de la fonction de transfert a partir de 
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la formule analytique, puis a en deduire une reponse 
impulsionnelle finie par transf ormee de Fourier inverse. 
On applique ensuite une operation de convolution pour le 
f iltrage . 

5 

Ainsi, en introduisant cette pre -compensation du champ 
proche a I'encodage, on definit une representation 
ambisonique modif iee (figure 8) , en adoptant comme 
representation transmissible des signaux exprim<§s dans le 
10 domaine frequentiel, sous la forme : 

Comme indique ci-avant, R est une distance de reference a 
15 laquelle est associe un effet de champ proche compense et 
c est la vitesse du son (typiquement 34 0 m/s dans l'air) . 
Cette representation ambisonique modifiee possede les 
memes proprietes de scalabilite (schematiqpiement 
representee par des donnees transmises "entour^es" pres de 
20 la fleche TR de la figure 1) et obeit aux memes 
transformations de rotation du champ (module 4 de la 
figure 1) que la representation ambisonique habituelle. 

On indique ci-apres les operations a mettre en oeuvre pour 
25 le decodage des signaux ambisoniques regus . 

On indique tout d 1 abord que 1 ' operation de decodage est 
adaptable a un dispositif de restitution quelconque, de 
rayon R 2 , different de la distance de reference R ci- 
30 avant, A cet effet, on applique des filtres de type 
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H^ FC( ^ /c ^ R/c) (co) , tels que decrits plus haut, mais avec des 
parametres de distance R et R 2 , au lieu de p et R. En 
particulier, il est a noter que seul le parametre R/c est 
a memoriser (et/ou transmettre) entre 1 1 encodage et le 
decodage . 

En se referant a la figure 12, le module de filtrage qui y 
est represente" ' est *'prevu par exemple dans une unite de 
traitement d'un dispositif de restitution. Les composantes 
ambisoniques regues ont ete pre-compensees a I 1 encodage 
pour une distance de reference R x en tant que seconde 
distance. Toutefois, le dispositif de restitution comporte 
une pluralite de haut s- par leurs disposes & une troisieme 
distance R 2 d f un point de perception auditive P, cette 
troisieme distance R 2 etant differente de la seconde 
distance precitee R x . Le module de filtrage de la figure 

12, sous la forme (co) , adapte alors, a la 

reception des donnees, la pre- compensation a la distance 
Ri pour une restitution a la distance R 2 . Bien entendu, 
comme indique ci-avant, le dispositif de restitution 
regoit aussi le parametre Ri/c. 

II est a noter que l r invention permet en outre de mixer 
plusieurs representations ambisoniques de champs sonores 
(sources reelles et/ou virtuelles) , dont les distances de 
reference R sont diff^rentes (le cas echeant avec des 
distances de reference inf inies et . correspondant a des 
sources lointaines) . Pref erentiellement , on filtrera une 
pre -compensation de toutes ces sources a une distance de 
reference la plus petite, avant de melanger les signaux 
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ambisoniques, ce qui permet a la restitution d'obtenir une 
definition correcte du relief sonore . 

Dans le cadre d ! un traitement dit de " focal i sat ion sonore" 
5 avec, a la restitution, un effet d 1 enrichissement sonore 
pour une direction choisie de 1 ' espace (a la maniere d'un 
projecteur lumineux eclairant dans une direction choisie 
en optique) , impliquant un traitement matriciel de 
focalisation sonore (avec pond^ration des composantes 
10 ambisoniques), on applique avantageusement 1 1 encodage de 
distance avec pre -compensation de champ proche de fagon 
combinee au traitement de focalisation. 

Dans ce qui suit, on d£crit un proc^de de decodage 
15 ambisonique, avec compensation du champ proche des hauts- 
parleurs, 3. la restitution. 

Pour reconstruire un champ acoustique encode suivant le 

formalisme ambisonique, a partir des composantes B% }1 et 
20 en utilisant des hauts -parleurs d'un dispositif de 
restitution qui prevoit un emplacement "ideal" d'un 
auditeur qui correspond au point de restitution P de la 
figure 7, 1 ' onde emise par chaque haut-parleur est definie 
par un traitement prealable de " re - encodage " du champ 
25 ambisonique au centre du dispositif de restitution, comtne 
suit . 

Dans ce contexte de "re-encodagre", on considere dans un 
premier temps et pour simplification que les sources 
30 emettent en champ lointain. 



41 

En se refer ant a nouveau a la figure 7, 1 1 onde emise par 
un haut-parleur d'indice i et d 1 incidence (0i et 8i) est 
alimente par un signal Si. Ce haut-parleur participe a la 

reconstruction de la composante B mn , par sa contribution 

Le vecteur c L des coefficients d'encodage associes aux 
hauts-parleurs d'indice i s 1 exprime par la relation : 



10 



Ci = 



Y 00 (&i>$i) 



[Bl] 



Le vecteur S des signaux emanant de 1 1 ensemble des N 
liauts-parleurs est donne par l 1 expression : 



15 



Si 



[B2] 
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La matrice d'encodage de ces N hauts-parleurs (qui 
correspond finalement a une matrice de "re-encodagre") , 
s 1 exprime par la relation : 



C= [c x C 2 ...C N ] 



[B3] 
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ou chaque terme c A represente un vecteur selon la relation 
[Bl] ci-avant. 

Ainsi, la reconstruction du champ ambisonique B 1 est 
definie par la relation : 



B = 



B 



II 



B 



= C.S 



[B4] 



La relation [B4] definit ainsi une operation de re- 
10 encodage, prealable a la restitution. Finalement, le 
decodage, en tant que tel, cons is te a comparer les signaux 
ambisoniques originaux et regus par le dispositif de 
restitution, sous la forme : 



15 



5 = 



B 00 



Bl 



[B5] 



aux signaux re -encodes B , pour def inir la relation 
genera le : 

B' = B [B6] 



20 
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II s'agit, en particulier, de determiner les coefficients 
d ! *une matrice de decodage D, qui verifie la relation : 

S = D.B [B7] 

De preference, le notnbre de hauts-parleurs est superieur 
ou egal au nombre de composantes ambisoniques a decoder et 
la matrice de decodage D s'exprime, en fonction de la 
matrice de re-encodage C, sous la forme : 

D = C T .(CC T y l [B8] 

ou la notation C T correspond a la transpos§e de la matrice 
C. 

II est §l noter que la definition d'un decodage verifiant 
des criteres differents par bandes de frequences est 
possible, ce qui permet d'offrir une restitution optimisee* 
en fonction des conditions d'ecoute, notamment pour ce qui 
concerne la contrainte de positionnement au centre 0 de la 
sphere de la figure 3, lors de la restitution. A cet 
effet, on prevoit avantageu semen t un filtrage simple, en 
egalisation f requentielle par paliers, §l chaque composante 
ambi sonique . 

Toutefois, pour obtenir une reconstruction d f une onde 
originellement encodee, il faut corriger I'hypothese de 
champ lointain pour les .hauts-parleurs, c ! est-a-dire 
exprimer 1' effet de leur champ proche dans la matrice de 
re-encodage C ci-avant et inverser ce nouveau systeme pour 
definir le d§codeur. A cet effet, en supposant une 
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concentricite des hauts-parleurs (disposes a une meme 
distance R du point P de la figure 7) , tous les hauts- 
parleurs ont un meme effet de champ proche F f \ J R/c) Uo) t sur 

chaque composante ambisonique du type B x % n . En 
introduisant les termes de champ proche sous la forme 
d'une matrice diagonale, la relation [B4] ci-avant 
devient : 

B>=Diag([\ F* /c (co) FX' c (a>) -F* /c (a;) F* /c (a>) -]).C.S [B9] 



La relation [B7] ci-avant devient 



r r 



S = D.Diag 



L l 



F R ' C (03) 



m 



(CO) 



3 



[BIO] 



Ainsi, l 1 operation de matrigage est prec^dee par une 
operation de filtrage qui compense le champ proche sur 

chaque composante B? nn , et qui peut etre mise en ceuvre 
sous forme numerique, comme decrit ci-avant, en reference 
a la relation [A14] . 

On retiendra qu'en pratique, la matrice C de "re-encodage" 
est propre au dispositif de restitution. Ses coefficients 
peuvent etre determines initialement par parametrage et 
caracterisation sonore du dispositif de restitution 
reagissant a un excitation predeterminee . La matrice de 
decodage D est, elle aussi, propre au dispositif de 
restitution. Ses coefficients peuvent etre determines par 
la relation [B8] . En reprenant la notation precedente ou 




45 



B est la matrice des composantes ambisoniques pre- 
compensees, ces dernieres peuvent etre transmises au 

dispositif de restitution sous forme matricielle B avec : 



B = Diag 



? Rlc 



RIc, 



3 



10 



Le dispositif de restitution decode ensuite les donn£es 
regues sous forme matricielle B (vecteur colonne des 
composantes transmises) en appliquant la matrice de 
decodage D aux composantes ambisoniques pre-compensees , 
pour former les signaux S A destines a alimenter les haut- 
parleurs HPi, avec : 



15 



20 



25 



S = 



f 5 0 

\ S NJ 



= D.B 



[Bii] 



En se referant a nouveau a la figure 12, si une operation 
de decodage doit §tre adaptee a un dispositif de 
restitution de rayon R 2 different de la distance de 
reference R x , un module d' adaptation prealable au decodage 
proprement dit et decrit ci-avant permet de filtrer chaque 

composante ambisonique B% n , pour 1 » adapter a un 
dispositif de restitution de rayon R 2 . L' operation de 
decodage proprement dite est effectuee ensuite, comme 
decrit ci-avant, en r£f6rence a la relation [Bll] . 
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On d§crit ci-apres une application de 1 ! invention a la 
synthese binaurale. 

On se refere & la figure 13A a sur laquelle un auditeur 
disposant d'un casque a deux ecouteurs d f un dispositif de 
synthese binaurale est represents. Les deux oreilles de 
1' auditeur sont disposees a des points respectifs O l 
(oreille gauche) et 0 R (oreille droite) de l'espace. Le 
centre de la t§te de 1' auditeur est dispose au point O et 
le rayon de la tSte de l 1 auditeur est de valeur a. Une 
source sonore doit etre pergue auditivement a un point M 
de l'espace, situe a une distance r du centre de la tete 
de 1 1 auditeur (et respectivement i. des distances r R de 
l l oreille droite et r L de 1 ! oreille gauche) . Par ailleurs, 
la direction de la source placee au point M est definie 

par les vecteurs r , r R et 7 L . 

De fagon g^n^rale, la synthase binaurale se definit comme 
suit- 

Chaque auditeur a xme forme d 1 oreille qui lui est propre . 
La perception d*un son dans l'espace par cet auditeur se 
fait par apprentissage, depuis la naissance, en fonction 
de la forme des oreilles (notamment la forme des pavilions 
et les dimensions de la tete) propre a cet auditeur. La 
perception d'un son dans l f espace se manifeste entre 
autres par le fait que le son parvient a une oreille, 
avant l 1 autre oreille, ce qui se traduit par un retard x 
entre les signaux a emettre par chaque ecouteur du 
dispositif de restitution appliquant la synthese 
binaurale . 
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Le dispositif de restitution est parametre initialement , 
pour un meme auditeur, en balayant une source sonore 
autour de sa tete, a une meme distance R du centre de sa 
5 tete. On comprendra ainsi que cette distance R peut §tre 
consideree comme une distance entre un "point de 
restitution" comme enonce ci-avant et un point de 
perception auditive (ici le centre O de la tete de 
1 1 auditeur) . 

10 

Dans ce qui suit, l'indice L est associg au signal a 
restituer par I'ecouteur accole a 1' oreille gauche et. 
l f indice R est associe au signal a restituer par ' t 
I'ecouteur accole k 1' oreille droite. En se referant a la -> 

15 figure 13B, on applique au signal initial S un retard pour : v| 
chaque voie destin^e a produire un signal pour un ecouteur 
distinct. Ces retards Tl, et x R sont f one t ion d'un retard 
maximum Xmax qui correspond ici au rapport a/c ou a, comme « 
indique precedemment , correspond au rayon de la tete de 

20 1' auditeur et c a la vitesse du son. En particulier, ces 
retards sont definis en fonction de la difference de 
distance du point O (centre de la tete) au point M 
(position de la source dont le son est a restituer, sur la 
figure 13A) et de chaque oreille a ce point M. 

25 Avantageusement , on applique en outre des gains respe.ctifs 
9Tl et g R , & chaque voie,. qui sont fonction d ! un rapport 
des distances du point O au point M et de chaque oreille 
au point M. Des modules respectifs appliques a chaque voie 
2 L et 2 R encodent les signaux de chaque voie, dans une 

30 representation ambisonique, avec pre -compensation de champ 
proche NFC (pour "Wear Field Compensation") au sens de la 
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presente invention. On comprendra ainsi que, par la mise 
en oeuvre du procede au sens de la presente invention, on 
peut definir les signaux issus de la source M, non 
seulement par leur direction (angles azimutaux G L et 8 R et 
5 angles d' elevation 8 L et 8 R ) , mais aussi en fonction de la 
distance separant chaque oreille r L et r R de la source M. 
Les signaux ainsi encodes sont transmis au dispositif de 

■restitution comportant .. des modules de . decodage 

ambisonique, pour chaque voie, 5 L et 5 R . Ainsi, on 
10 applique un encodage/d^codage ambisonique, avec 
compensation de champ proche, pour chaque voie (ecouteur 
gauche, ecouteur droit) dans la restitution avec synthese 
binaurale (ici de type "B-FORMAT") , sous forme dedoublee. 
La compensation de champ proche s'effectue, pour chaque 
15 voie, avec comme premiere distance p une distance r L et r R 
entre chaque oreille et la position M de la source sonore 
a restituer. 



On decrit ci-apres une application de la compensation au 
20 sens de 1' invention, au contexte de 1 • acquisition sonore 
en representation ambisonique. 

On se refere a la figure 14 sur laquelle un microphone 141 
comporte une pluralite de capsules transductrices, 

25 capables de capter des pressions acoustiques et restituer 
des signaux electriques S x , ...,S N . Les capsules CAPi sont 
agencees sur une sphere de rayon r predetermine (ici, une 
sphere rigide, telle qu'une balle de ping-pong par 
exemple) . Les capsules sont espacees d'un pas regulier sur 

30 la sphere. En pratique, on choisit le nombre N de capsules 
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en fonction de l'ordre M desire pour la representation 
ambisonique . 

On indique ci-apres, dans le contexte d'un microphone 
comportant des capsules agencees sur une sphere rigide, 
comment compenser 1 1 ef f et de champ proche, des l ! encodage 
dans le contexte ambisonique. On montrera ainsi que la 
pre- compensation du champ proche peut s'appliquer non 
seulement pour la simulation de source virtuelle, comme 
indique ci-avant, mais aussi a 1 1 acquisition et, de fagon 
plus generale, en corabinant la pre- compensation de champ 
proche & tous types de traitements impliquant une 
representation ambisonique. 

En presence d'une sphere rigide (susceptible d'introduire 
une diffraction des ondes sonores regues) , la relation 
[Al] donnee ci-avant devient : 



00 

p r&0= £ ~ 72737, : Z<« y m«fe) [Cl] 
w=0 (Ar) h m (kr) 0£n ^ m 

cr=±l 



Les derivees des fonctions de Hankel spheriques h" m 
obeissent a la loi de recurrence : 



(2m + Y)h m (x) = rn (*)-(«+!) h m +l( x ) [C2] 



On d£duit les coraposantes ambisoniques B? nn du champ 
initial a partir du champ de pression a la surface de la 



10 



15 



25 



# 
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sphere, en mettant en oeuvre des operations de projection 
et d' egalisation donnees par la relation : 

BZ„=EQ m <p r \Y£ n > 4rr IC3] 



Dans cette expression, EQ m est un filtre egaliseur qui 
compense vine ponderation W m qui est liee a la directivite 
des capsules et qui inclut en outre la diffraction par la 
sphere rigide. 

L 1 expression de ce filtre EQ m est donnee par la relation 
suivante : 



Les coefficients cie ce filtre d 1 egalisation ne sont pas 
stables et on obtient un gain infini en tres basses 
frequences. D'ailleurs, il convient de noter que les 
composantes harmoniques spheriques, elles-m§mes, ne sont 
20 pas d l amplitude finie lorsque le champ sonore n'est pas 
limite a une propagation d'ondes planes, c'est-a-dire 
issues de sources lointaines , comme on 1 x a vu 
precedemment . 



Par ailleurs, si, plutot que de pr^voir des capsules 

encas trees dans une sphere solide, on prevoit des capsules 

de type cardioxdes, avec une directivite en champ lointain 
donnee par 1 1 expression : 
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G(6) = a + (1 - a) cos 6 [ C5 ] 



En considerant ces capsules montees sur un support 
"transparent acoustiquement" , le terme de ponderation a 
5 compenser devient : 

W m =j m {ajm{kr)-j{l-a)jm\krj) [C6] 

II apparait encore que les coefficients d'un filtre 
10 d 1 egalisation correspondant a 1 1 inverse analytique de 
cette ponderation donnee par la relation [C6] sont 
divergents pour les trSs basses frequences. 

De fagon generale, on indique que pour tout type de 
15 directivite de capteurs, le gain du filtre EQ m pour 
compenser la ponderation W m liee & la directivity des 
capteurs est infini pour les basses frequences sonores. En 
se r^ferant a la figure 14, on applique avan t ageus emen t 
une pre -compensation de champ proche dans l r expression 
20 mSme du filtre d' egalisation EQ m , donnee par la relation : 



BQ^'*W-*%&j*L [C7] 

Ainsi, les signaux 3 X a S N sont r^cuperes du microphone 
141. Le cas echeant, on applique une pre-egalisation de 
ces signaux par un module de traitement 142. Le module 143 
permet d'exprimer ces signaux dans le contexte 
ambisonique, sous forme matricielle. Le module 144 
applique le filtre de la relation [C7] aux composantes 
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ambisoniques exprimees en fonction du rayon r de la sphere 
du microphone 141. La compensation de champ proche 
s'effectue pour une distance de reference R en tant que 
seconde distance. Les signaux encodes et ainsi filtres par 
5 le module 144 peuvent etre transmis, le cas echeant , avec 
le parametre representatif de la distance de reference 
R/c. 

Ainsi, il apparait dans les differents modes de 

10 realisation lies respect ivement a la creation d ! une source 
virtuelle en champ proche, a 1 1 acquisition de signaux 
sonores issues de sources reelles, ou meme & la 
restitution (pour compenser un effet de champ proche des 
hauts-parleurs) , que la compensation de champ proche au 

15 sens de la presente invention peut s'appliquer a tous 
types de traitements faisant intervenir une representation 
ambisonique. Cette compensation de champ proche permet 
d'appliquer la representation ambisonique a une 
multiplicity de contextes sonores oxi la direction d'une 

20 source et avantageusement sa distance doivent etre prises 
. en compte. De plus, la possibility de la representation de 
phenomdnes sonores de tous types (champs proches ou 
lointains) dans le contexte ambisonique est as sure e par 
cette pre -compensation, du fait de la limitation a des 

25 valeurs reelles finies des composantes ambisoniques . 

Bien entendu, la presente invention ne se limite pas a la 
forme de realisation decrite ci-avant a titre d r exemple ; 
elle s'etend a d'autres variantes . 
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Ainsi, on comprendra que la pr6- compensation de champ 
proche peut £tre int<§gree, a l'encodage, autant pour une 
source proche que pour une source lointaine. Dans ce 
dernier cas (source lointaine et reception d'ondes 
5 planes) , la distance p exprimee ci-avant sera consideree 
comme infinie, sans modifier de fa<?on substantielle 
1' expression des filtres H m donnee ci-avant. Ainsi, le 
traitement utilisant des processeurs d'effet de salle qui 
fournissent en general des signaux decorreles utilisables 
10 pour modeliser le champ diffus tardif (reverberation 
tardive) peut etre combine a une pre -compensation de champ 
proche- On peut considerer que ces signaux sont de meme 
energie et correspondent a une part de champ diffus;^. 

correspondant a la composante omnidirective ^ — Bqq^ 
15 (figure 4) . On peut alors construire les diverses *f 
composantes harmoniques spheriques (avec un ordre M' 5" 
choisi) en appliquant une correction de gain pour chaque "t 
composante ambisonique et on applique une compensation de" 
champ proche des hauts-parleurs (avec une distance de 
20 reference R separant les haut-parleurs du point de 
perception auditive comme represents sur la figure 7) . 

Bien entendu, le principe d'encodage au sens de la 
presente invention est generalisable a des modeles de 
25 rayonnement autres que des sources monopolaires (reelles 
ou virtuelles) et/ou des hauts-parleurs. En effet, toute 
forme de rayonnement (notamment une source etalee dans 
l'espace) peut etre exprimee par integration d ! une 
distribution continue de sources 61ementaires ponctuelles . 
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En outre, dans le contexte de la restitution, il est 
possible d' adapter la compensation de champ proche a tout 
contexte de restitution. A cet effet, il peut etre previa 
de calculer des fonctions de transfert (re-encodage des 
composantes harmoniques spheriques de champ proche pour 
chaque haut-parleur , compte tenu d'une propagation- reelle 
dans la salle ou le son est restitue) , ainsi qu'une 
inversion dece re-encodage pour redefinir le decodage. 

On a decrit ci-avant un procede de decodage dans lequel on 
appliquait un systeme matriciel faisant intervenir les 
composantes ambisoniques . Dans une variante, il peut etre 
prevu un traitement generalise par trans formees de Fourier 
rapides (circulaire ou spherique) pour limiter les temps 
de calcul et les ressources informatiqu.es (en terme de 
memoire) necessaires au traitement de decodage. 

Comme indique ci-avant en reference aux figures 9 et 10, 
on constate que le choix d'une distance de reference R par 
rapport a la distance p de la source en champ proche 
introduit une difference de gain pour differentes valeurs 
de la frequence sonore. On indique que le proced<§ 
d'encodage avec pre -compensation peut etre couple a une 
compression audionumerique permettant de quantifier et 
d'ajuster le gain pour chaque sous-bande f requentielle . 

Avantageusement, la presente invention s f applique a tous 
types de systemes de spatialisation sonore, notamment pour 
des applications de type "realite virtuelle" (navigation 
dans des scenes virtuelles dans l'espace tridimensionnel, 
jeux avec spatialisation sonore tridimensionnelle, 
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conversations de type "chat" sonorisees sur le reseau 
Internet), a des sonif ications d 1 interfaces, a des 
logiciels d 1 edition audio pour enregistrer, mixer et 
restituer de la musique, mais aussi a 1 1 acquisition, §l 
5 partir d' usage de microphones tridimensionjiels, pour la 
prise de son musicale ou cinematographique, ou encore pour 
la transmission d 1 ambiance sonore sur Internet, par 
exemple pour des "WebCam" _ sonorisees . 
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REVEHD I CATIONS 

1. Precede de traitement de donnees sonores, dans lequel : 

a) on code des signaux representatif s d'au moins un son se 
propageant dans 1 1 espace tridimensional et issu d f une 
source situee a une premiere distance (p) d'un point de 
reference (O) , pour obtenir une representation du son par 
des composantes - (B^) -exprimees .dans .une . base 
d ' harmoniques spher iques , d 1 origine correspondant audit 
point de reference (O) , 

b) et on applique auxdites composantes (B^) une 
compensation d'un effet de champ proche par un filtrage 
qui est fonction d'une seconde distance (R) definissant 
sensiblement, pour une restitution du son par un 
dispositif de restitution, une distance entre un point de 
restitution (HPi) et un point (P) de perception auditive. 



2. Proced6 selon la revendication 1, dans lequel, ladite 
source etant lointaine du point de reference (O) , 

- on obtient des composantes d'ordres successifs m pour la 
representation du son dans ladite base d 1 harmoniques 
spheriques, et 

- on applique un filtre (1/FJ dont les coefficients, 
appliques chacun a une composante d'ordre m, s 1 expriment 
analytiquement sous la forme de 1 ' inverse d'un polynome de 
puissance m, dont la variable est inversement 
proportionnelle k la frequence sonore et a ladite seconde 
distance (R) , pour compenser un effet de champ proche au 
niveau du dispositif de restitution. 
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3. Procede selon la revendication 1, dans lequel, ladite 
source etant une source virtuelle prevue a ladite premiere 
distance (p) , 

- on obtient des composantes d'ordres successifs m pour la 
representation du son dans ladite base d'harmoniques 
spheriques, et 

- on applique un filtre global (H m ) dont les coefficients, 
appliques chacun £ une compos ante d'ordre m, s ' expriment 
analytiquement sous la forme d'une fraction, dont : 

- le numerateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la 
frequence sonore et a ladite premiere distance (p) , 
pour simuler un effet de champ proche de la source 
virtuelle, et 

- le denominateur est un polynome de puissance m, 
dont la variable est inversement proportionnelle 3. la 
frequence sonore et a ladite seconde distance (R) , 
pour compenser 1 ! effet du champ proche de la source 
virtuelle dans les basses frequences sonores . 

4. Proced(§ selon l'une des revendications precedentes, 
dans lequel on transmet au dispositif de restitution les 
donnees cod<§es et f iltrees aux etapes a) et b) avec un 
parametre representatif de ladite seconde distance (R/c) . 

5. Procede selon l'une des revendications 1 a 3, dans 
lequel, le dispositif de restitution coraportant des moyens 
de lecture d' un support memoire, on memorise sur un 
support memoire destine a etre lu par le dispositif de 
restitution les donnees codees et filtr€es aux etapes a) 
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et b) avec un parametre representatif de ladite seconde 
distance (R/c) . 

6. Procede selon 1'une des revendications 4 et 5, dans 
lequel, prealablement a une restitution sonore par un 
dispositif de restitution comportant une pluralite de 
hauts-parleurs disposes a une troisieme distance (R 2 ) 
dudit point de perception auditive (P) , on applique aux 
donnees codees et £iltr€es un filtre d' adaptation 
(H ra (R1/C/R2/c) ) dpnt les coefficients sont fonction desdites 
seconde (R x ) et troisieme distances (R 2 ) - 

7. Procede selon la revendication 6, dans lequel les 
coefficients dudit filtre d' adaptation (Hm <Ri/c,R2/o } , 
appliques chacun a une composante d'ordre m, s'expriment 
analytiquement sous la forme d'une fraction, dont : 

- le numerateur est un polyn6me de puissance m, dont la 
variable est inversement proportionnelle & la frequence 
sonore et a ladite seconde distance (R) , 

- et le denominateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la frequence 
sonore et a ladite troisidme distance (R 2 ) • 

8. Procede selon l'une des revendications 2, 3 et 1, dans 
lequel, pour la mise en oeuvre de I'etape b) , on prevoit : 

- pour des composantes d ! ordre m pair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 
d'ordre deux ; et 

- pour des composantes d'ordre m impair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 
d'ordre deux et une cellule supplementaire d'ordre un. 
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9 . Procede selon la revendication 8 , dans lequel les 
coefficients d'un filtre audionumerique, pour une 
composante d ! ordre m, sont definis a partir des valeurs 

5 numeriques des racines desdits polynomes de puissance m. 

10. Procede selon 1'une des revendications 2, 3, 1, 8 et 
9, dans lequel lesdits polyn6mes sont des polyn6mes de 
Bessel . 

10 

11. Procede selon I'une des revendications 1, 2 et 4 d. 10, 
dans lequel on prevoit un microphone comport ant un reseau 

• de transducteurs acoustiques agenc^s sensiblement sur la. 
1 surface d'une sphere dont le centre correspond^ 
15 sensiblement audit point de reference (O) , pour obtenir^ 
lesdits signaux representatif s d f au moins un son se^ 
propageant dans 1 1 espace tridimensionnel . £ 

12. Proc£d6 selon la revendication 11, dans lequel on 
20 applique a l'Stape* b) un filtre gl9bal pour, d'une part, 

compenser un effet de champ proche en fonction de ladite 
seconde distance (R) et, d' autre part, 6galiser les 
signaux issus des transducteurs pour compenser une 
ponderation de directivite desdits transducteurs. 



25 



13. Procede selon l'une des revendications 11 et 12, dans 
lequel on prevoit un notnbre de transducteurs fonction d'un 
nombre total choisi de composantes pour representer le son 
dans ladite base d' harmoniques spheriques . 
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14. Precede selon l'une des revendications precedences, 
dans lequel on choisit a l'etape a) un nombre total de 
composantes dans la base des harmoniques spheriques pour 
obtenir, & la restitution, une region de I'espace autour 

5 du point de perception (P) dans laquelle la restitution du 
son est fidele et dont les dimensions sont croissantes 
avec le nonibre total de composantes . 

15. Procede selon la revendication 14, dans lequel on 
10 prevoit un dispositif de restitution comportant un nombre 

de haut-parleurs au moins 6gal audit nombre total de 
composantes . 

16. Procede selon l'une des revendications 1 a 5 et 8 a 
15 13, dans lequel : 

- on prevoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une 
distance choisie d ! un auditeur, 

- on obtient, pour cet auditeur, une information de 
20 ressenti de la position dans 1 1 espace de sources sonores 

situees a une distance de reference predetermine (R) de 
1 1 auditeur , et 

- on applique la compensation de l'etape b) avec ladite 
distance de reference sensiblement en tant que seconde 

25 distance . 

17. Procede selon l'une des revendications 1 a 3 et 8 a 
13, prises en combinaison avec l'une des revendications 4 
et 5, dans lequel : 
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- on prevoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une . 
distance choisie d f un auditeur, 

on obtient, pour cet auditeur, une information de 
5 ressenti de la position dans l'espace de sources sonores 
situees a line distance de reference predeterminee (R 2 ) de 
1* auditeur, et 

- prdalablement a une restitution^ sonore par le dispositif 
de restitution, on applique aux donnees codees et fil trees 

10 aux etapes a) et b) un filtre d' adaptation (H m (R ^ ,R2/c) ) 
dont les coefficients sont fonction de la seconde distance 
(R) et sensiblement de la distance de reference (R 2 ) • 

18. Proced<§ selori l'une des revendications 16 et 17, dans . 

•ft 

15 lequel : ^ 

- le dispositif de restitution comporte un casque k deux ^ 
ecouteurs pour les oreilles respectives de l 1 auditeur, et 

- separement pour chaque Scouteur, on applique le codage 
et le filtrage des etapes a) et b) pour des signaux 

20 respect if s destines a alimenter chaque ecouteur, avec, en 
tant que premiere distance (p) , respect ivement une 
distance (r^r^) separant chaque oreille d'une position 
(M) d'une source a restituer. 

25 19. Procede selon l'une des revendications precedentes, 
dans lequel on met en forme, aux Stapes a) et b) , un 
systeme matriciel comportant au moins : 

- vine matrice (B) comportant lesdites composantes dans la 
base des harmoniques spheriques, et 
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- une matrice diagonale (Diag(l/F m )) dont les coefficients 
correspondent a des coefficients de filtrage de l'etape 
b) , 

et on multiplie lesdites matrices pour obtenir une matrice 
5 resultat de composantes compensees {B). 

20. Procede selon la revendication 19, dans lequel : 

- Le dispositif de restitution comporte une pluralite de 
haut-parleurs disposes sensiblement a une meme distance 
(R) du point de perception auditive (P) , et 

- pour decoder lesdites donnies codees et filtr^es aux 
e tapes a) et b) et former des signaux adapt es pour 
alimenter lesdits haut-parleurs : 

* on forme un systeme matriciel comportant ladite 

. matrice resultat (2?) et une matrice de d€codage (D) 
pr^determinee , propre au dispositif de restitution, 
et 

* on obtient Tine matrice (S) comportant des 
coefficients representatif s des signaux 
d 1 alimentation des hauts-parleurs par multiplication 

de la matrice des composantes compensees { B ) par 
ladite matrice de decodage (D) . 

21. Dispositif d 1 acquisition sonore, comportant un 
25 microphone muni d'un reseau de transducteurs acoustiques 
disposes sensiblement sur la surface d'une sphere, 
caracterise en ce qu'il comporte en outre une unit6 de 
traitement agencee pour : 

- recevoir des signaux 6manant chacun d'un transducteur , 
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- appliquer auxdits signaux un codage pour obtenir une 
representation du son par des composantes (Bm^) exprimees*. 
dans line base d 1 harmoniques spheriques, d'origine 
correspondant au centre de ladite sphere (O) , 

- et appliquer auxdites composantes (B^ 0 ) un filtrage qui 
est fonction, d'une part, d r une distance correspondant au 
rayon de la sphere (r) et, d 1 autre part, d'une distance de 
reference -(-R) ■- - 

22. Dispositif selon la revendication 21, caractgrise en 
ce que ledit filtrage consiste, d'une part, a 4galiser, en 
fonction du rayon de la sphere, les signaux issus des 
transducteurs pour compenser une ponderation de 
directivite desdits transducteurs et, d 1 autre part, a 
compenser un effet de champ proche en fonction d'une 
distance de reference choisie (R) , def inissant 
sensiblement, pour une restitution du son, une distance 
entre un point de restitution (HP ± ) et un point (P) de, 
perception auditive. 
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